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We theoretically evaluate changes in the magnetic potential arising from the magnetic field near 
superconducting thin films. An example of an atom chip based on a three-wire configuration has 
been simulated in the superconducting and the normal conducting state. Inhomogeneous current 
densities within the superconducting wires were calculated using an energy-minimization routine 
based on the London theory. The Meissner effect causes changes to both trap position and oscillation 
frequencies at short distances from the superconducting surface. Superconducting wires produce 
much shallower micro traps than normal conducting wires. The results presented in this paper 
demonstrate the importance of taking the Meissner effect into account when designing and carrying 
out experiments on magnetically trapped neutral atoms near superconducting surfaces. 
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I. INTRODUCTION 

Atom chips with microfabricated and nanofabricated 
field-generating elements are useful devices for the co- 
herent manipulation of ultracold atoms A diversity 
of potentials with high spatial resolution can be gener- 
ated near the chip surface using state-of-the-art fabri- 
cation technology of metals, semiconductors, supercon- 
ductors and optical waveguides. By means of these po- 
tentials Bosc-Einstein condensates and degenerate Fermi 
gases are routinely prepared. Significant progress in co- 
herent manipulation and state-sensitive detection of sin- 
gle atoms has been achieved in recent years d, Q • Ap- 
plications of atom chips include atom interferometers 
HjSHQi ultra sensitive magnetic field sensors [1,3 and 
portable experimental systems for quantum gases [10[ . 

A key role on atom chips is played by the atom-surface 
interactions, such as undesirable spin-decoherence mech- 
anisms and attractive Casimir-Polder forces [H, [3 ■ 
The lifetime of magnetic traps near the chip surface 
is limited by decoherence mechanisms produced by the 
near-field noise radiation from thermally induced cur- 
rents in conductive surfaces [ll|, [13, [13, [IB • Cooling the 
chip can reduce those thermal currents and thus increase 
the lifetime of the magnetic traps. An important increase 
in the lifetime of several orders of magnitude is expected 
when the surface layer crosses the transition from the nor- 
mal to the superconducting state [l^. The achievement 
of such conditions promises a new class of experiments 
with cold atoms integrated with nanostructured surfaces 
that will allow coherent control over the atoms into the 
submicron range. Exciting proposals for coupling atoms 
to nanodevices such as mechanical oscillators |17| or su- 
perconducting circuits [isj outline new perspectives on 
experimental research at the interface between atomic 
and solid-state physics. 

Although the usage of superconducting microstruc- 
tures on atoms chips was already proposed in 1995 by 



Wcinstcin and Libbrech t Il9l , and later on in several the- 
oretical studies 0, [13, l2l| |. the impact of the Meissner 
effect [22I, [2^ on the magnetic trapping potential has not 
been considered as yet. In this paper, we show that the 
magnetic-field expulsion from the superconductor has an 
important effect on the magnetic trap. Our theoretical 
results will be essential for the proper design of super- 
conducting atom chips, a subject on which the first ex- 
perimental results have recently been reported [l^, [IB] . 

Accurate methods to calculate the potential landscape 
near the chip surface is a prerequisite for cold-atom ex- 
periments on atom chips. In the case of using supercon- 
ducting wires, these calculations must take into account 
the Meissner effect, which implies that the intensity of in- 
duced or applied currents decays exponentially into the 
interior of the superconductor with a penetration depth 
A, which is also the depth to which the magnetic field 
penetrates the superconductor. Magnetic-field calcula- 
tions are especially important at short distances from 
the chip surface, where atoms are not easily accessible 
by imaging methods and where the Meissner effect can 
have an important impact on the trap parameters. 

The main goal of this paper is to investigate how the 
Meissner effect modifies magnetic traps near supercon- 
ducting thin films. To accomplish this task, magnetic 
fields and trap parameters such as position, oscillation 
frequencies and trap depth were calculated in simulated 
chips containing thin-film wires. Simulations were car- 
ried out for the vortex-free superconducting state and for 
the normal conducting state, and the differences between 
the two cases arc analyzed. 

Current distributions in superconducting thin films 
were calculated in the frame of the London theory [2^ . 
Despite its fundamental nature, exact solutions of the 
London theory exist only for trivial cases such as a single 
sphere or a single cylinder in a homogeneous magnetic 
field [12 . Numerical methods are therefore necessary to 
calculate current density distributions in thin-film mi- 
crostructures. Brandt and Mikitik [26| reported on how 
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lus increases linearly in the radial directions: 




FIG. 1: Representation of a typical atom-chip geometry that 
provides three-dimensional magnetic confinement [2^. This 
theoretical example is used in this paper to study the prop- 
erties of magnetic micro traps near superconducting thin 
films. Three parallel wires on the chip surface generate a 
two-dimensional quadrupole field B2D that provides radial 
confinement. The width of the quadrupole wires and the sep- 
aration between them are denoted by w and v, respectively. 
Underneath the chip surface there are two offset wires in the 
perpendicular direction to supply longitudinal confinement. 
The offset wires are located sd, z = q/2 and z = — g/2. 



to obtain numerical solutions of the London theory for 
strips with rectangular cross section in a perpendicular 
homogeneous magnetic field and/or with applied electric 
current. More general geometries can be solved using 
commercial programs which, however, have severe lim- 
itations. For example, most of them provide accurate 
solutions only if the thickness h of the thin film is similar 
to the penetration depth A [27| . In the present paper we 
overcome this limitation and provide an algorithm that 
provides accurate solutions of the London theory by find- 
ing the current distribution that minimizes the free en- 
ergy. A similar minimization method [2^ has been used 
to obtain the magnetization curves of arrays of super- 
conducting strips in homogeneous magnetic fields. The 
numerical method presented in this paper can solve more 
general geometries in arbitrary inhomogeneous magnetic 
fields, including most of the geometries that are typically 
present in atom chips. 



II. MAGNETIC CONFINEMENT ON AN ATOM 

CHIP 

A magnetic micro trap can be realized with the atom- 
chip geometry represented in Fig. [1] Three parallel thin- 
film wires on the chip surface generate a two-dimensional 
confining field JB2D (magnetic guide). The current Ic in 
the central wire is opposite in direction to the currents 
Ibi and Ib2 in the two outer wires. The magnetic guide 
forms at the position (a;o,yo)j where the magnetic field 
of the central current Ic is canceled by the bias field 
generated by Ibi and Ib2- The field B2D forms a two- 
dimensional quadrupole field around (a;o,J/o)- Its modu- 



\B2d\ = a\/{x - xq)^ + iy- VoY 



(1) 



Here, a is the gradient of the magnetic guide in the ra- 
dial directions. For simplicity, these three parallel wires, 
which we will refer to as quadrupole wires, are assumed 
to be identical in size. The width and the thickness of 
each quadrupole wire are denoted by w and h, respec- 
tively. The outer wires are separated from the central 
wire by v. 

Longitudinal confinement is achieved by means of the 
inhomogeneous offset field Bq generated by two offset 
wires perpendicular to the quadrupole wires. The two 
offset wires are located below the chip surface, and sep- 
arated from the quadrupole wires by k. They are driven 
with identical currents /q. 

The offset field Bq is superimposed onto the two- 
dimensional confining field B2d, hi such a way that a 
magnetic trap forms between the two offset wires around 
the point (xq, i/q, 0). Near the centre of the trap. 
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where oq and are the fist- and second-partial deriva- 
tives of Bo with respect to the corresponding directions, 
evaluated at (a;o,yo,0). The trap forms only if bz > 0; 
thus for distances yo+k smaller than about 0.6 q. The off- 
set field bo at the trap center normally suffices to prevent 
Majorana spin-flip transitions [l^l and the consequent 
loss of atoms in Bose-Einstein condensates. However, in 
experiments with thermal clouds at temperatures of sev- 
eral micro K, it might be necessary to externally apply an 
additional homogeneous offset field Bo,ext = (0, 0, bo,ext)- 
The centre of the trap 



ao{bo + bo^ext) „ 
xo,yo T— — 5 — '0 



(3) 



is slightly displaced from the position (xo^yoyO) of the 
magnetic guide. 

The offset field changes the radial potential from lin- 
ear to parabolic, with a harmonic oscillation parameter 
a'^ /{bo + bo,ext)- The magnetic trap is then characterized 
by the radial and longitudinal oscillation frequencies 



m{bo + bo.ext 
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Here, gp is the Lande factor, /is is the Bohr magneton, 
mp is the magnetic quantum number, and m is the atom 
mass. 

The y component of Bo causes a small rotation of the 
longitudinal axis of the magnetic trap. In the particular 
case that the magnetic trap is in the plane a; = 0, and 
thus Ibi = Ib2, the quadrupole field can be expressed as 
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Fig. [2|). The current density is assumed homogeneous 
within each thin strip, although it may vary from strip 
to strip. The free energy of this system is the sum of 
the magnetostatic energy of the currents and the kinetic 
energy of the electrons [1^, and can be written in the 
form 



3Af -iN 



(6) 



FIG. 2: Sketch of the division used to calculate the applied 
current distributions in the quadrupole wires. A cross section 
of the central wire is shown. Every wire is divided up into a 
large number of longitudinal strips with squared cross section 
of side s. The current density within each strip is assumed 
homogeneous. 



and the rotation occurs about the y axis |2£ 
angle 9 = ao/a. 



with 



III. CALCULATION OF INHOMOGENEOUS 
CURRENT DENSITIES IN 
SUPERCONDUCTING THIN FILMS 

Current densities are homogeneous in normally con- 
ducting wires. Inhomogeneous current densities in the 
superconducting wires were calculated numerically using 
the London theory. General solutions of the London the- 
ory for the microstructure shown in Fig. [1] can be ob- 
tained by linear superposition of two separate cases. The 
first of these describes the behavior of the microstructure 
when the three quadrupole wires carry the respective cur- 
rents Ic, Ibi and Ib2- The second case describes the 
behavior of the microstructure when each offset wire is 
driven with a current Iq and no current is applied to the 
quadrupole wires. In the second case, induced screening 
currents in the quadrupole wires can modify the param- 
eters of a magnetic trap formed near the chip surface. 

For simplicity, the offset wires are assumed to be in- 
finitely thin. This approximation is valid if the width of 
the offset wires is much smaller than fc, in which case 
neither the screening currents in the quadrupole wires 
nor the magnetic fields near the chip surface depends on 
the current distribution in the offset wires. The homo- 
geneous offset field Bo.ext = (0, 0, feo.ert) is not distorted 
by the Meissncr effect because the longitudinal demagne- 
tizing factor of a strip quickly tends to zero as the strip 
length increases to infinity [31| . 



where /„ is the electric current along the strip n and 
Mnrn IS the mutual inductance between the strips n and 
TO. The general expression for Mnm is js^ . Issj 
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where Snm is the Kronecker delta, J„ = In/ is the 
current density, and r„ denotes the position of point 
within the strip n. The first and second terms repre- 
sent the magnetic and the kinetic inductances, respec- 
tively. The integrals are carried out over the volumes 
of the corresponding strips. Since the current density is 
homogeneous within each strip, the magnetic term can 
be approximated by classical formulas tabulated in Ref. 
[3^ . The integral of the kinetic term has a trivial solu- 
tion. The matrix elements then become 
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where / is the length of the wires and d is the distance 
between the centers of the considered strips. 

The superconducting current density is obtained by 
finding the set of values ..,3Af that minimize the 

fimction F{Ii, I2, Isn)- This is accomplished with the 
method of the Lagrange multipliers, and imposing that 
the total current flowing in each wire is flxed. The con- 
straints are 
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The equations to find the superconducting currents are 
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A. Applied currents in the quadrupole wires 
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iV+l,...,27V,(10) 



Current density distributions in the superconduct- 
ing quadrupole wires are calculated with an energy- 
minimization numerical procedure. Every quadrupole 
wire is divided up into a large number N of thin lon- 
gitudinal strips with squared cross section of side s (see 
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where Ai , A2 and A3 are the Lagrange multipliers. The 
solution of this system of linear equations is a set of values 
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{In}n=i,...,3N that represent the current distribution in 
the superconducting wires. 

For long wires (/ > lOw) the calculated current distri- 
bution does not depend on I. This limit is valid in all 
the examples shown in this paper, where all wires are 
assumed infinitely long. 

Low values of s improve the accuracy of the solution 
at the expense of long computation time. We have found 
that all values lower than A/2 produce practically the 
same numerical results, and therefore, s = A/2 is in gen- 
eral a good choice for calculations. 

Calculating the mutual inductance Mmn between two 
strips is not incompatible with the general definition of 
mutual inductance for two closed circuits. One can imag- 
ine that every quadrupole wire is part of a closed circuit 
that includes the current drivers and the wires between 
the chip and the drivers. The mutual inductance between 
two strips can be defined as the contribution of the two 
strips to the total mutual inductance between the closed 
circuits of which they form part. 
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B. Screening currents in the quadrupole wires 

Screening currents arise in superconductors in the pres- 
ence of external magnetic fields as a consequence of 
the Mcissncr effect. Screening currents can be calcu- 
lated using the energy-minimization method described 
in this section. This numerical method requires to de- 
compose the screening currents in small current elements 
{In}n=i....,N- In a magnetostatic situation, the screen- 
ing currents are closed, and therefore, a decomposition 
based on small magnetic dipoles or small closed current 
elements is the most appropriate. 

The geometry of the closed current elements /„ de- 
scribed in this paragraph is suitable to evaluate the en- 
ergy and the flux including the kinetic term [s^l. The 
superconducting body is divided up into small cubes of 
side s. The current is assumed to be homogeneous within 
every small cube, although both the intensity and the 
direction of the current might vary from cube to cube. 
Closed current elements /„ similar to magnetic dipoles 
are built by grouping the cubes in sets of four, in the 
way illustrated in Fig. [3l The centers of the four cubes 
lie in the same plane. Sets of four cubes can be built 
in planes parallel to the x axis {x sets), to the y axis {y 
sets) or to the z axis {z sets). In this manner, any cube 
that is not on the wire surface belongs to twelve different 
sets: four x sets, four y sets and four z sets. Sets are 
sorted with numbers. Every set has associated a current 
/„, n being the number of the set. The current element 
/„ is distributed within the set n in a way that is similar 
to a magnetic dipole, as shown in Fig. [31 The current 
/„ changes its direction by 90 degrees as it passes from 
a cube to the next cube of the set. The direction of /„ 
in each cube is described by the unit vectors C/„_i, C/„_2, 
Un,3 and C/„,4. 

The mutual inductance Mnm between two sets n and 



FIG. 3: Sketch of the sets of cubes in which the quadrupole 
wires are divided to calculate the screening currents. The 
lower part illustrates how the current 7„ is distributed within 
the set n. The direction of the current, which is different 
in every cube, is represented by four unit vectors Un,i, Un,2 
Un;i and Un,i- The effective surface of the set is symbolized 
by a gray square of side y/2s. 



m is obtained by summing the contributions of the mu- 
tual inductances between the separate cubes of both sets. 
Assuming one set n made up of the cubes ni, 712 , and 
?i4 and another set m made up of the cubes mi, TO2, TO3 
and 7714, the total mutual inductance between the two 
sets is 



Um,b) : 



(11) 



where 
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is the contribution of the cubes Ua and mb to Mnm- 
The scalar product Un.a • Um,b accounts for the fact that 
this contribution depends on the angle between the cur- 
rent directions. Here, r„.a denotes the position of point 
within the cube ria, and J„ = In/ is the current density 
in the set n. The integrals are carried out over the vol- 
ume of the corresponding cubes. The double i nteg ral of 
the magnetic term was calculated numerically |34l |. and 
the integral of the kinetic term has a trivial solution. The 
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matrix elements are then approximated by 
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where d is the distance between cube centers. 

For the reasons mentioned above, the mutual induc- 
tance between two cubes is a senseless physical idea un- 
less each cube is regarded as part of a close circuit. To 
understand the meaning of Mnm, the screening current 
tubes can be regarded as a collection of closed circuits 
with a certain inductance matrix. In this way the mu- 
tual inductance between two cubes can be understood as 
the contribution made by the cubes to the total mutual 
inductance between the two current tubes in which the 
cubes are included. This idea also applies to the self- 
inductance. 

Every set of cubes is also characterized by the effective 
surface Sn, which is represented in Fig. [3] by the gray 
area. The effective surface is defined so that the product 
Mnmlm IS the flux produccd by through Sn- The 
modulus of Sn is 2s^, and its direction is determined by 
the right-hand rule from the direction of the current /„. 
Following the notation of Fig. [31 the effective surface can 
be expressed by Sn = 2s^C/„_i x Un,2- 

The solution of the London theory for this system is 
the current distribution that cancels the flux -including 
both the magnetic and the kinetic terms- in the interior 
of the superconducting wires. It is possible to demon- 
strate [22| that this condition is equivalent to the free- 
energy minimization performed in the previous section. 
The equations to be solved are formulated so that the 
flux through every set of cubes is null: 



N 



^ ^ MnmJm + Sn ' -B„ — . 



1,...,7V. (14) 



Here, N is the total number of sets, Mnm symbolizes 
the elements of the inductance matrix, B„ is the external 
magnetic field at the position of the set n, and Sn is 
the effective surface of the set n. The first term in this 
equation represents the total flux <f>5,„ induced by the 
screening currents onto the set n. The second term is 
the magnetic flux $0,71 of the external field onto the set 
n. 

Due to its large size, the matrix Mnm cannot be in- 
verted using any of the mathematic programmes that 
are generally available. As an alternative, an itera- 
tive method has been used to solve Eq. (fT4|) . In the 
fist step, a homogeneous screening current distribution 

is assumed: Im = -Sm ■ Bm/Mmm] m = l,...,N. 
This distribution will not satisfy Eq. (fTi)) . and the 
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Sm=i -^'^nmlrn crcatcd by thc assumed 



flux <I> 

screening currents onto each set n will not cancel thc 
flux created by the offset wires $o.n- In the second 



(2) 

step, the current distribution is calculated by /A = 
In^ - in\'^n,n + $5;|j)/M„„, which wiU generate a flux 
that will be more similar to the desired solution. The pro- 
cess continues until the convergence condition is satisfied: 
'J'o.n + <f>5 |j < 10~^$o,n, after e iterations. The number 
^ is a convergence factor that has no physical meaning. 
This value is chosen by trial and error. The best choice 
depends on the geometry of the superconducting body. 
If ^ is too high, the method is not convergent; but if ^ is 
too low, the convergence is very slow. ^ can vary from 
set to set and also from step to step. The particular val- 
ues of ^ used for our calculations will not be shown here 
since thc results do not depend on them and since there 
arc many other possibilities. 

Once the closed current elements ..,7V have 

been obtained, the total current in every single cube is 
calculated by summing the contributions of the sets of 
which the cube forms part. From the total currents in 
the single cubes, the magnetic field outside the supercon- 
ducting body can be calculated. 

Thc method presented in this section can be used 
to calculate screening currents in superconducting thin- 
films induced by arbitrary inhomogeneous magnetic 
fields. In the example shown in Fig. (TJ the magnetic field 
generated by the offset wires has no component along the 
X direction. For that reason, only y sets and z sets are 
used. 



IV. NUMERICAL RESULTS 

In this section we assess how the Meissner effect al- 
ters the magnetic-trap parameters in thc superconduct- 
ing microstructure depicted in Fig. [1] In all the examples 
shown in this paper, the width of the quadrupole wires 
and the separation between them are w = 5 /.im and 
w = 5 ^m, respectively. The penetration depth A in the 
superconducting wires is 100 nm, which is a typical value 
for metallic superconductors. The offset wires are at a 
distance of fc = 5 /im underneath the quadrupole wires. 
The thickness of thc quadrupole wires h as well as the 
separation between the two offset wires q will be varied 
in the examples in order to demonstrate their effect on 
the magnetic-trap parameters. 



A. Quadrupole magnetic guide generated by the 
superconducting chip 

First wc analyze thc position of thc magnetic guide in 
the X, y plane. Thc magnetic guide can be positioned 
within a large area above thc chip surface by changing 
the ratios 



Ic 



Ibi 



'B2 



/3 = 7^. (15) 



Figure m illustrates the trajectories corresponding to con- 
stant a (dashed lines) and to constant (3 (solid lines) 
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FIG. 4: Position {xo,yo) of the magnetic guide in the x,y- 
plane as a function of the ratios a and (3 for the supercon- 
ducting and the normal conducting chip. Solid Hues are tra- 
jectories generated by varying a while keeping /3 constant. 
Dashed lines are trajectories generated by varying /3 while 
keeping a constant. The values of a and /3 are written on 
the corresponding trajectory. The width and the thickness of 
the quadrupole wires, and the separation between them are 
w = 5/^m, h = and v — 5/im, respectively. 
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FIG. 5: Radial gradient obtained for different positions of 
the magnetic guide {xo,yo) keeping the sum of the currents 
Is = Ic + Ibi + Ib2 at a constant value of 1 A. a) Radial 
gradient ajvc in the normal conducting chip, b) ratios of the 
gradient asc in the superconducting chip to the gradient aMC 
in the normal conducting chip. The width and the thickness 
of the quadrupole wires, and the separation between them are 
w = 5jj,m, h = 150/im and v = 5fj,m, respectively. 

for the superconducting and the normal conducting chip. 
Differences between the two cases are noticeable when the 
distance between the magnetic trap and the chip surface 
is smaller than the width of the wires. 

In principle, both the position in the x,y plane and 
the radial gradient of the quadrupole field B2D depend 




FIG. 6: Current density Jc in the central wire for two different 
values of the thickness h. The section of the wire is shown. 
The width of the wires and the separation between them are 
w = 5/im and v — 5/im, respectively. The applied currents 
are Ic=65fj.A, and Jsi=/s2=467/iA (a=0.07; 13=1). 



on the applied currents Ic, Ibi and Ib2- Once the ra- 
tios a and (3 have been chosen to position the magnetic 
guide, the radial gradient can be varied by changing the 
value of Is = Ic + Ibi + Ib2- Fig. El^a) shows for con- 
stant =1 A the radial gradient qnc in the normal 
conducting chip as a function of the position of the mag- 
netic guide. The radial gradient for other values of Is can 
be obtained by linear scaling. For the superconducting 
case, the gradient asc was calculated in the same way, 
keeping Is at a constant value of 1 A. Figure [5l[b) shows 
the ratio asc/o-NC- Superconducting wires produce con- 
siderably lower radial gradients than normal conducting 
wires. The radial gradient of B2D is related with the 
radial oscillation frequency of the micro trap by Eq. 

Changes in the trapping field caused by the Meiss- 
ner effect become more pronounced when the supercon- 
ducting wires are thicker or when the penetration depth 
is smaller. Either thinner wires or longer penetration 
depths imply more homogeneity in the superconducting 
current densities, which results in magnetic fields which 
are more similar to those produced by normal conductors. 
Figure [6] shows the current-density distribution Jc{x,y) 
along the central wire for two different thicknesses. Three 
regimes can be distinguished. If /i ^ A, the current den- 
sity decays exponentially from the surface and shows a 
sharp peak in each corner. If ft, A, the current density 
becomes homogeneous along the y axis, having two max- 
ima at X = w/2 and x = —w/2. For extremely thin wires, 
the kinetic energy gets so high that the current density 
becomes almost homogeneous, allowing the magnetic fiux 
to penetrate the film. 

In the case of normal conducting wires, the magnetic- 
trap parameters were independent of the thickness h. 
This is illustrated by comparing the numerical results ob- 
tained for different values of h. The variations in the x, y 
position and in the radial gradient produced by varying 
h between 50/im and SOO/^m were, respectively, less than 
0.01/im and less than 0.1% at any position within the 
area represented in Figs. |4]and[5l On the contrary, our 
numerical calculations demonstrate that the magnetic- 
trap parameters depend considerably upon the value of 
the thickness h when the chip is superconducting. For 
example, while for /i=150 nm the radial gradient in the 
superconducting chip is 0.4 times the radial gradient in 
the normal conducting chip at 0.5 fim from the chip sur- 



z{nm) 

FIG. 7: Results obtained for the magnetic field Bo when 
the offset wires are driven with equal currents /o = lmA and 
no current is applied to the quadrupole wires, a) stream 
lines of the induced screening currents integrated along the 
y-direction in the central wire. The plotted current density is 
20fiA/\ine. b) Field lines in the plane a; = above the super- 
conducting chip, c) Field lines in the plane x — above the 
normal conducting chip. The field lines indicate the direction 
of the field but the density of lines does not show the intensity 
of the field. One can appreciate the expulsion of the magnetic 
field from the interior of the superconducting wires. Calcula- 
tions are performed for the following geometrical parameters: 
w = 5/^m, V = 5^m, /i=150nm, k — 5fim and q = 30/^m. 
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FIG. 8: Longitudinal component of the magnetic field along 
the z-direction at {xo,yo) = (0,2/^m). Three different cases 
are represented: normal conductor, superconductor with h — 
150nm, and superconductor with h = 350nm. Jo^lmA. The 
other geometrical parameters are the same than in Fig. [T] 

face above the central wire (see Fig. [5]), this reduction 
factor is 0.6 for /ii=50 nm and 0.3 for /i2=500nm at the 
same position. Therefore, the thickness of the thin-film 
wires becomes relevant in the superconducting state. 



B. Longitudinal confinement in the 
superconducting state 

The analysis presented in this section is restricted to 
magnetic traps located in the plane a; = 0. Different 
distances from the chip surface will be considered. 

Figure [7] shows the screening currents in the central 
quadrupole wire as well as the magnetic field lines of 
the offset field Bq for the superconducting and the nor- 
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FIG. 9: This figure compares the trap parameters bo and 
bz between the superconducting and the normal conducting 
chip, a) Ratio bo,sc/bo,NC- b) Ratio bz,sc/bz,NC- The hor- 
izontal axis, in logarithmic scale, represents the distance q 
between the two offset wires, and the vertical axis, in linear 
scale, represents the position of the magnetic trap yo- Data 
are represented in the plane a;=0. The other geometrical pa- 
rameters are: w = 5/^m, i) = 5fim, h = 150nm, k = bfim. 
The region in which no trap forms is left of the dashed line. 



mal conducting chip. Figure [8] represents the longitu- 
dinal component of the magnetic field calculated along 
the z direction at {xo,yo) = (0,2 ^m) for three differ- 
ent cases: normal conductor, superconductor with h = 
150 nm, and superconductor with h = 350 nm. As with 
the results obtained for B2d, differences between the su- 
perconducting and the normal conducting states become 
larger with increasing h. The screening currents in the 
superconducting quadrupole wires reduce the z compo- 
nent of the magnetic field Bq at the positions z = —9/2 
and z = q/2. This effect entails a decrease of the trap 
depth along the longitudinal direction. This reduction 
is of about 15% at 2 fim from the surface, and becomes 
higher than 25% at distances of 1 fim or shorter. For 
above than 10 ^m, the reduction is lower than 5%. 

The parameters bo, gq and bz that describe the inho- 
mogcneous offset field Bq were numerically calculated for 
the superconducting chip (SC) and the normal conduct- 
ing chip (NC) as a function of q and t/q. Figure [DJa) 
shows the ratio bo^sc/bo,NC- The horizontal axis repre- 
sents the distance q between the two offset wires, and the 
vertical axis represents the position of the magnetic trap 
yo- As observed in this figure, the Meissner effect in the 
superconducting wires slightly increases the value of 60 ■ 
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FIG. 10: Trap parameter ao in the normal conducting chip 
(a) and in the superconducting chip (b). 7o=lmA. The hor- 
izontal axis, in logarithmic scale, represents the distance q 
between the two offset wires, and the vertical axis, in linear 
scale, represents the position of the magnetic trap j/o- Data 
are represented in the plane x=0. The other geometrical pa- 
rameters are: w = 5fim, v = 5/^m, h = 150nm, k = 5fim. 
The region in which no trap forms is left of the dashed line. 

This increase becomes more significant as the distance 
between the two offset wires q gets longer, and the mag- 
netic trap gets closer to the surface. Figure [9jb) shows 
the ratio hz,sc /i>z,NC- The longitudinal oscillation fre- 
quency of the microtrap is related to hz by means of Eq. 

As seen in Fig. [Hlthe longitudinal frequency is dra- 
matically reduced by the Meissner effect when the offset 
wires are close to each other. For high values of q the ef- 
fect is the opposite, and the longitudinal frequencies are 
slightly higher in the superconducting chip. 

Figure [TUl compares the value of oq between the super- 
conducting and the normal conducting chips. The pa- 
rameter oq is related with the angle of rotation of the trap 
as explained in Sec. [Ill The calculated values of oq were 
significantly lower in the superconducting microstructure 
than in the normal conducting microstructure. 



V. SIMULATION OF A MAGNETIC MICRO 
TRAP 

In this last section we apply the numerical results pre- 
sented in Sec. IIVI to a typical example of a magnetic 
micro trap. Figure [11] shows the isopotential curves of 
a magnetic trap generated by the atom chip depicted in 
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FIG. 11: Isopotential curves of a magnetic trap generated 
by the atom chip shown in Fig. [T] in the superconducting 
(SC) and in the normal conducting (NC) state. The ap- 
plied currents are the same in both cases: Ic ~ 0.2mA, 
Ibi = Ib2 = 1.4mA (q = 0.0714, /3 = 1 and Jg = 3mA) 
and lo = 2mA. Homogeneous offset field bo,ext = 25/iT is 
externally applied to stabilize the micro trap against Majo- 
rana losses. Following the notation of Fig. [T] the geometry 
of the microstructure is described hy w = 5^m, v = 5/im, 
h = 150nm, q = lOO/xm and k — 5/im. The penetration 
depth is lOOnm in the superconducting wires. The two upper 
graphs show the isopotential curves in the plane y — yo- The 
dashed lines in the two upper graphs represent the longitudi- 
nal axis of the micro trap, which is rotated about the y-axis as 
explained before. The two lower graphs show the isopotential 
curves in the plane perpendicular to the chip surface along the 
longitudinal axis. The dashed lines in the two lower graphs 
represent the chip surface. The parameters of this micro trap 
are represented in Table 1. The magnetic field changes by 
4/iT per contour. 

Fig. [T] in the superconducting and in the normal con- 
ducting state. The applied currents are the same in both 
cases. The relevant trap parameters are summarized in 
Table [Tj The micro trap forms closer to the surface in 
the superconducting chip than in the normal conduct- 
ing chip. In the present example, the Meissner effect 
produces an important reduction in the radial oscillation 
frequencies as well as a slight increase of the longitudinal 
oscillation frequencies, as predicted by Figs. [5|and[9l 

The most remarkable feature of the supercoducting 
chip is a significant decrease of the trap depth towards 
the surface, which is a consequence of the reduction of a 
shown in Fig. [5] In the shown example, the trap depth 
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is reduced by about 80% in the superconducting chip. 





sc 


NC 


a (Tm~") 


6.9 


8.4 


ao (Tm-i) 


0.2 


0.3 


60 (mT) 


2.5 


2.3 


b, (mTm-2) 


5700 


5240 


yo {^J■n-^) 


2.0 


2.3 


wi (s-i) 


2-K- 95 


2-K- 92 




2-K- 1650 


27r- 2020 


e 


1.7° 


2.1° 



TABLE I: Parameters of the micro trap shown in Fig. Illl for 
the superconducting (SC) and the normal conducting (NC) 
states. Oscillation frequencies have been calculated for *^Rb. 



come relevant when the distance between the micro trap 
and the superconducting surface is smaUer than the 
width of the wires. The most dramatic consequence of 
the Meissner effect is a significant reduction of the trap 
depth. In general, the Meissner effect has to be taken 
into account when designing and carrying out experi- 
ments with neutral atoms magnetically trapped near su- 
perconducting surfaces. Although the results shown in 
this paper have been obtained for the specific example il- 
lustrated in Fig. [U these conclusions can be generalized 
to any atom chip made with superconducting thin films. 
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It has been numerically solved in Mathematica using the 
function NIntegrate. 

The energy and flux can be properly evaluated from the 
closed current elements {/n}n=i,...,jv when these occupy 
the whole volume of the superconducting body and when 
the overlapping between neighbor current elements is to- 



tal, in the sense that a homogeneous distribution of closed 
current elements generate null electric current at any in- 
ternal point that is not on the surface of the supercon- 
ducting body. 



